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ABSTRACT 


Wc report on the gamma-ray observations of giant molecular clouds Orion A and B 
with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. 
The gamma-ray omission in the energy band between ~ 100 MeV and ~ 100 GeV is 
predicted to trace the gas mass distribution in the clouds through nuclear interactions 
between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production 
ernes- section for the nuclear interaction is known to ~ 10% precision which makes the 
LAT a powerful tool to measure the gas mass column density distribution of molecular 
clouds for a known CR intensity. We present here such distributions for Orion A and B, 
and correlate them with those of the velocity integrated CO intensity (Woo) at a 1° x 1° 
pixel level. The correlation is found to be linear over a Wco range of ~ 10 fold when 
divided in 3 regions , suggesting penetrAtiou of nuclear CRs to most of the cloud volumes. 
The Wco-to-mass conversion factor, A'co, is found to be ~ 2.3x 10 20 cm -2 (K km s“ l )" 1 
for the high-longitude part of Orion A (I > 212°), ~ 1.7 times higher than ~ 1.3x 10 20 
found for the rest of Orion A and B. We interpret the apparent high Xco ' n the high- 
longitude region of Orion A in the light of recent works proposing a non-linear relation 
between Hj and CO densities in the diffuse molecular gas. Wqo decreases faster than 
the Hj column density in the region making the gas “darker’ to Wco- 

Subject headings molecular clouds: general molecular clouds: individual(Orion A. 
Orion B) 


1. Introduction 

The Orion A and B clouds are the archetypes of local giant molecular clouds (GMCs) where 
interstellar gas condenses and stars are formed (e.g., Bergin k Tafalla 2007; Bally 2008, and ref- 
erences therein). The clouds have been studied in various wavebands including millimeter obser- 
vations of the transition lines between CO rotational states, especially from J = 1 to J — 0 (e.g., 
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Sanders et al. 1984; Maddalcna et al. 1986; Dame et al. 1987, 2001; Wilson ct al. 2005; Fukui et al 
2011), infrared emission (e g., Beichman et al. 1988), attenuation of star light (e.g., Dobashi et al. 

2005), and near infrared extinction (Row lee k Froebrich 2009; Froebrich k Rowlca 2010; Dobashi 
2011). The two clouds are prime targets for the Large Area Telescope (LAT) on-board the Fermi 
Gamma-ray Space Telescope (Fermi) in the research of molecular clouds and CR interaction be- 
cause they lie isolated from. the Galactic plane and do intense gamma ray point source overlaps 
with the clouds (Abdo et al. 2009d; Abdo et aL 2010b). 

Gamma rays from the Orion Monoceros region were first detected by COS-B in the energy 
range between 100 MeV and 5 GeV (Caraveo et al. 1980; Bloemen et al. 1984). EGRET detected 
gamma rays in the range between 100 MeV and ~ 10 GeV (Digel et al. 1995, 1999). In these 
studies, the gamma- ray intensity distribution in a region including Orion A, B and Monoceros R2 
was fitted with three independent contributions, one proportional to the atomic hydrogen (H I) 
column density, another proportional to the CO line intensity (Wco) 1 . and the last, a presumed 
isotropic distribution. Under the assumptions that Wco traces the Hj column density, the CR 
spectrum doesn’t change in the region and H T spin temperature (7s) is constant, the ratio Xco 
was determined^, from the ratio of the gamma-ray intensities associated with the H I and CO 
distributions, to be Xco — (2.6 ± 1.2) x 10 20 cm 2 (K km a -1 ) -1 (Blocmen et al. 1984) and Xco 
(1.35 ± 0.15) x 10 20 cm 2 (K km 8 J ) -1 (Digcl et al. 1999). The ratio was not separately measured 
for the three clouds, Orion A, B and Monoceros R2, due to the limited statistics and spatial 
resolutipn of the instruments. We note that Strong et al. (1988) determined Xco on the diffuse 
Galactic gamma rays observed by COS-B to be Xco = (2.3 ±0.3) x 10 20 cm' 2 (K km s“ 1 )"' and 
Dame et al. (2001), by comparing smoothed infrared intensity and Wco distributions across the 
Galaxy, determined it to be Xco = (1-8 ±0.3) x 10 20 cm" 2 (K km a" 1 ) -1 . 

Since the publications on the EGRET data (Digel et al. 199o ; 1999), much progress has been 
made in studies on Orion A and B: new observational data became available (e.g., Dame et al. 2001; 
Lombardi k Alves 2001; Wilson et al. 2005; Kalberla et al. 2005; Dobashi et al. 2005; Rowlos k Froebrich 
2009; Froebrich k Rowlos 2010; Dobashi 2011); study of the molecular clouds was renewed (e.g., 
Wilson et al. 2005; Bally 2008); a new modeling of the Galactic diffuse gammarray emission was pro- 
posed incorporating large-scale CR propagation (Strong k Moskalenko 1998; Strong ot al. 2000); 
theoretical calculations of collisional CO rotational-level excitation were revisited (Mengcl et al. 


1 We define Wco as the velocity- integrated intensity of the transition line between J = 1 to J = 0 in 

3 Our Xco a facto* converting VVco to mam column dunnity measured in unite of the proton mass in cloud 
concentrations predominantly consuming of H a .' In some literature Xoo * osed as tlm factor converting Wco to 
Ilj column density. Where Woo traces Hj accurately and the chemical state of hydrogen is predominantly in Ha, 
the 2 definitions are expected to agree. Tho helium and heavier atoms are assumed to be mixed uniformly in the 
Interstellar gas with the solar abundance. Wc warn reader*) that comparison of Xoo vahwe calculated on different 
CO surveys and gamma-ray observations are not straightforward due to differences in their calibration procedure 
(e.g see Bronfman et al. 1988, for the CO calibration factor) as well as in the assumptions on the CR composition 
and the associated croae sections. 



2001; Flower 2001; Cecchi-Pcstellini et aJ. 2002; Balakrishnan et al. 2002; Wernli et al. 2006; Sheplcr ct id. 
2007; see also Kalberla et al. 2005; Liszt 2006, 2007) and the distance to the Orion nebula in the 
Orion A cloud was measured accurately (Sandstroin et al. 2007; Menton ct al. 2007; Hirota et al. 

2007; Kim et al. 2008). 

The Fermi Gamma-ray Space Telescope mission, launched on 2008 June 11, has been sur- 
veying the sky with the Large Area Telescope (LAT) since 2008 August. Its wide field of view, 
large effective area, improved spatial resolution, and broad energy coverage provide much higher 
sensitivity relative to its predecessor EGRET (Atwood et al. 2009; Abdo et al. 2009a). 

Studies based on EGRET observations have established that gamma rays from Galactic molec- 
ular clouds are dominated by neutral piou decays (which wr refer to as the ^pionic gamma rays” or 
“pionic emission" ) in the energy band bet ween 0.2 GeV and 10 GeV (Bertsch nt al. 1993; Digel et al. 

1995, 1999). Orion A and B are located far (~ 8.8 kpc) from the Galactic center 3 and displaced 
from the Galactic plane by — 140 pc. The two clouds are only ~ 400 pc away from the solar 
system where spectra of CR species upto the sub-TeV domain are predicted to be similar to those 
measured directly at the Earth after correction for the solar modulation. 

We can uow analyze Orion A and B through the high-energy gamma rays detected by the Fermi 
LAT in the light of the recent, developments and study the relation between IFco ^d mass column 
density (or Xco) in various parts of the Galaxy and obtain the total mass of the clouds 4 . The 
improved spatial resolution and higher gamma-ray statistics provided by the Fermi -LAT allow us 
to determine the relation on angular scales of 1 x 1 deg 2 (pixels), without being directly aifccted by 
the thermodynamical, chemical, or radiation environment insido the Orion clouds, albeit within the 
limited angular resolution of t.h« Fermi LAT and uncertainties due to any unresolved weak sources 
and CR flux variation. The results can be used conversely to study various environmental effects 
on Xco in the translucent parts of clouds where most gas in Orion A and B resides and where the 
Xqo factor has not been straightforward to derive (e.g., van Dishocck Black 1986; Magnani et aL 
1988; Bolatto ct al. 1999; Magnani et al. 2003; Bell et al. 2006; Snow & McCall 2006, Bell et al. 

2007; Burgh et al. 2007; Wall 2007; Sheffcr et al. 2008). 

Theoretical analyses have long suggested that A'co depends on the environment and the Wco* 
A^H?) relation may be nonlinear (e.g., Kutncr & Leung 1985; Dickman ct al. 1986; Maloney L Black 
1988; Taylor et al. 1993; Bolatto et al. 1999; Magnani et al. '2003; Bell ct aL 2007; Burgh et al. 
2007). Suggestions have also been made that Xco depends on the relative abundances of CO, C I. 
and C II (e.g. van Dishoeck U Black 1988; Hollenbach ct al. 1991; Kopp et al. 2000). The existence 
of gas not. traced by H I and CO at the interface between the two phases (the “dark gas") has been 


'We aiamme thr distance between the Sun and the Galactic center to be 8 A kpc and the Galactic rotation velocity 
near the Sun to be 230 km a" 1 . 

*The miM of Orion A and B la distributed mostly in the column density range corresponding to a “translucent 
cloud whoac lina-of-sight visual attenuation (Ay) is typically between . and 5 mag and has u(Hj) typically between 
100 and 2000 cm -5 (e.g. van Disboedc k Black 1988). 




discovered (Grenier et al. 2005; Ade et «J. 2011). The relation between the fraction of carbon in CO 
and H 2 density in translucent and diffuse clouds has been updated based on observations and nu- 
merical simulations, for example, by Burgh et al. (2010); Wolfire «t al. (2010). Glover et al. (2010). 
Our results will be interpreted in the light of these recent works. The Wco W(Hj) relation will be 
characterized including the “dark gas.” and the measured mass column density will be related to 
the A\ value at which the relation is predicted to become non-linear. 

In this paper we analyze diffuse gamma rays spatially associated with the molecular clouds'’ 
Orion A and B, extract their pionic gamma-ray components, obtain mass distributions, and compare 
them with those predicted for Wco measured by Fukui et al. (2011) and Dame et al. (2001). In 
Section 2 we describe the gamma-ray event selection applied in this analysis. The analysis procedure 
is described in Section 3 in 4 subsections: the spatial templates used to extract mass column density 
associated with multiple emission components are given in Subsectiou 3.1; energy-binned spatial 
fits on the templates are described in Subsection 3.2; the pionic emission is extracted from the 
spectra obtained in the spatial fits and Xco Is calculated thereon in Subsection 3.3; and the total 
H 2 masses of Orion A and B are estimated in Subftection 3.4. In Section 4, we assess systematic 
uncertainties in the analyses: check the Xco results with recent infrared excess emission maps by 
Dobash: (2011); summarize the results; and interpret them in the light of recent studies of the 
relatiou between the H*j and CO fraction In the translucent clouds The paper is concluded in 
Section 5. 


2. Observations and Data 

The data used In this analysis were obtained in the nominal all-sky survey mode between 
2008 August 4 and 2010 March ll 6 . We select events classified as Pass 6 Diffuse class which has a 
high gamma- ray purity (Atwood et al. 2009). Among the events, wc limit the reconstructed zenith 
angle to be less than 105° to greatly reduce gamma rays coming from the limb of the Earth’s 
atmosphere. We select the good time intervals (GTIs) of the observations by excluding events that 
were taken while the instrument rocking angle was larger than 52°. Another cut is made on the 
reconstructed gamma-ray energy at E m , n - 178 MeV and * = 100 GcV to reduce systematic 
uncertainty of the LAT effective area and residual background events induced by CRs. Gamma 
rays in a rectangular region of 30° x 30° centered at (/ = 2L0°, 6 = -20°) are then selected for 
later analyses. We refer to the region as the regiou-of-intercat (ROI) and the set of events as the 
data set. 


‘By molecular cloud* we mean spatially identified clouds without distinguishing the small admixture of atomic 
and ionisad hydrogens therein. 

6 Mission Elapsed Time 230,557,413 b through 290.000.000 s where zero is art at 00 00 UTC on 2001 January 1. 
During the period, the LAT was operated in the survey mode with the rocking angle 35 deg (2008 August 4 to 2009 
July 9), 39 deg (2009 July 9 to 2009 September 3) and 50 deg (2009 September 3 to 2010 March 11). 




The data set consists of 1,132,436 events of which 901,929 are between 178 MeV and 1 GeV, 
224,753 between 1 GeV and 10 GeV, and 5,754 hetwoen 10 GeV and 100 GeV. They are binned in 
150 x 150 equal-area pixels ( Hammer- A i toff projection) in Galactic coordinates with 0.2° gridding 
an their reconstructed arrival directions, and in 22 logarithmic bins between = 178 MeV and 
£ ^, t = 100 GeV on their reconstructed energies 

The map of counts integrated over the energy range of the data set is shown in Fig. 1. We can 
visually identify Orion A and B near the center of the region and the outer Galactic plane in the 
upper part We note that Monoceros R2 is also visible between Orion A/B and the outer Galactic 
plane. 


3. Analyses 

The analyses presented here begin by finding the relationship between the spatial distributions 
of gamma rays and Wqq, the most widely used proxy of Ha, in the Orion clouds and by studying the 
proportionality between the two and its spatial dependence within the Orion clouds. The analyses 
proceed in 3 steps. 

In the first step, the spatial distribution of thp “background* gamma rays, i.e., the gamma 
rays not associated with the H 2 clouds, is determined by using spatial distribution templates, for 
the H 1 gas, for the inverse Compton scattering (IC) component, for the point sources, and for a 
presumed isotropic component (Subsections 3.1). We then fit, in Subsection 3.2, the gamma-ray 
spatial distribution in each of the 22 energy bins as a sum of the “background" distribution and a 
distribution tentatively associated with the H 2 gas (Hj- template). The “background" is subtracted 
from the measured gamma-ray intensity distribution and the remainder is defined as the gamma ray 
intensity distribution associated with the H 2 gas with which is correlated pixel-by-pixel. We 
note that the gamma-ray intensity measures the mass column density in the Hj gas for a known 
CR spectrum. We repeat the fit with 2 alternative Ha-templates. 

In the second step (Subsections 3.3 and 3.4), the energy-binned gumma-ray cmiaaivity for the 
H 2 gas (J3, in cq. (1)) are assembled as the gamma-ray spectrum for each of the 3 ^-templates. 
The spectrum is then fitted os a sum of the gamma rays produced in the piouic and bremsstrahlung 
processes. 

In the third step, the gammarray intensity distribution associated with the pionic emission is 
converted to the mass column density. The IVco-mass conversion factor ( Xqo ) i* calculated via 
two methods, one by comparing the gammar-ray counts associated with the H I gas and with the H 2 
gas (the H 2 /H I method) and the other by dividing the gamma-ray counts of the pionic emission 
by the number of pionic gamma rays expected per unit gas mass (the pionic method). In the 
first method, we assume the CR spectrum is uniform in the local H I region within Galactoccntric 
radius of 8 - 10 kpc (see Subsection 3.1,1) and in the Orion clouds. In the latter method, we 
assume the CR spectrum including its absolute flux is known in the Orion clouds. We validate 



these assumption using GALPROP. 

We use GALPROP (Strong & Moskalenko 1998; Strong et a!. 2000) with the parameter set 
labeled as GALDEF 54.77Xvarh7S. This parameter set is known to reproduce reasonably well the 
Galactic diffuse gamma-ray emission observed with the LAT (Abdo et al. 2009f) 7 . We refer to the 
results obtained by running GALPROP with this parameter sot as the GALPROP results in this 
paper. 


3.1. Spatial distribution templates 

Initially wo assume the gamma-ray emission from the ROI to be made of 4 “background" 
components and one “signal" Hj component . each emitting gamma rays with a characteristic spatial 
distribution. The 4 “background” components are spatially associated with the diffuse H ! gas, the 
inverse Compton (IC) scattering by electrons 8 off interstellar radiation fields, the point sources, and 
the sum of extragalartic diffuse emission (including unresolved sources) and backgrounds inducted 
by CRs in the instrument. We assume the last surn to be isotropic and refer to it as the isotropic 
component. We ignore the contribution from ionized hydrogen gas (H 11) because its density is low 
when averaged in 1 x 1 deg 2 pixels (< 0.5cm ~ 8 ) and its total mass is negligible in the ROI (Gordon 
1969; O’Dell 2001). 

All spatial components except for the IC component arc assumed to have, individually, an 
energy- independent, underlying spatial distribution in Galactic coordinates (1,6). Another impor- 
tant underlying assumption is that the nuclear CR spectrum is uniform over the ROI. We make 
spatial templates for the 22 energy bins by convolving the spatial distributions with the energy- 
dependent point spread function (PSF) and exposure for the individual energy bins. Hence the 
spatial templates are cuorgy dependent. In Subsection 3.3 we will show that the spectra of the 
gamma-ray emissions associated with the H I and H 2 gas consist of the pionic and brcmsstrahlung 
components. 

The gamma-ray intensity / 7 (/, 6) for the i-tb energy bin is interpreted as the sum of the five 
contributions, each being the product of the normalization factor for the i-th energy bin and the 
spatial template. 

Ml, 6) = A,N(H m.b) + B,/V(H,)(»,6) + /C,(1,6) + + A. (1) 

i 

The normalization factors are: A, for the H 1 gas; I3 t for the component associated with clouds 
consisting predominantly of H 2 ; IC{ for the inverse Compton component; C t} 6i } for the j-th point 


: A detailod description of GALDEF files can he found at hUp://gnlprop.*tanfonl.odu 
6 We refer to electrons as a sum of e 4 and e 



source at and /}< for the isotropic component which is assumed not to depend on (/, b). The 

normalization factors we determined independently for the 22 energy bins. We note that ICi are 
fixed at the values given by GALPItOP, because the Bpatial distribution is highly correlated with 
the isotropic component, and the IC component is sub-dominant in the ROI. 

Later in Subsection 3.2, wc will explore 3 templates for Ha, two based on Wco and one on 
Wco plus the “dark gas" proposed by Grcnior et al. (2005). 


3.1.1. Diffuse H / gas template 

Atomic hydrogen gas (H 1) is broadly distributed in the Galaxy with a total mass exceeding 
that of molecular hydrogen (Hj) (c.g., Fcrriirc 2001; Snow & McCall 2006). In the outer Galaxy 
where the Orion clouds are located, the mass column density of H I is lower than that of Ha at the 
Orion clouds (Kalberla et al. 2010, 2005) 

We used the Leiden /Argentine/ Bonn (LAB) survey data (Kalberla et al 2005) corrected for 
optical thickness by adopting a constant spin temperature (Ts) of 125 K as the H I gas spatial 
distribution template (sec Fig. 2a). The LAB intensity distribution is divided into five annuli cen- 
tered at the Galactic center as has been done in other Fermi diffuse emission analyses (Abdo et al. 
2010c). Their inner and outer Galactoccutric radii ( R ) are: 8 to 10, 10 to 11.5, 11.5 to 16.5, 16.5 
to 19, and 19 to 50 kpc. The line-of-sigh*. velocity distribution of the H 1 gas in the Orion region 
overlaps that of the CO gas associated with the Orion clouds and that of the local H I annulus 
(R = 8-10 kpc) quite well. 

Gamma-ray contributions from all the H I annuli overlapping our ROI have been included 
in the analyses. In the fitting, the CR intensity is treated independently at each annulus. The 
contributions from annuli other than the local one (R — 8 - 10 kpc) are through the periphery 
of the LAT PSF and less than ~ 5% in gamma-ray counts. Hence our analyses are insensitive to 
variation in the CR intensity and/or spectrum among the neighboring annuli. 

The spin temperature of H I gas, Ts, is not well constrained in the region nor known to be 
uniform over the ROI: its quoted value in the literature ranges between ~ 90 K and ~ 400 K 
(e.g. Mohan ct al. 2004a, b). We estimate, later iu this paper, the contribution to the overall 
systematic error from this uncertainty by repeating the analysis for Ts = 250 K and 90 K. No 
significant concentration of cold H 1 is known wound Orion A and B at lwge-scale (> 1 x 1 deg 2 ) 
(Kalberla et al. 2010). An exploratory study of cold H I mixed in selected H 2 cloud cores has found 
the mean H I fraction to bo less than 0.5% (Krto et al. 2008). So we can safely ignore such a 
mixture in the analysis. 

Gamma rays are produced in the H I gas through the pionic wid bremsstrahlung processes 
with intensities proportional to the CR nuclear and electron spectra in the gas, respectively. 



8.1.2. Molecular cloud template 


We try 3 H 2 templates to represent the H 2 spatial distribution in the ROI. In making the 
templates we assume that the H 2 column density is proportional to Woo measured by two CO ( J — 
1 — > 0) surveys, one from NANTEN (Fukui et ai. 2011) covering the areas around the Orion clouds 
with effective resolution of 4' and the other being the Galactic survey by Dame et al. (2001) covering 
the ROI with angular resolution of 8.7*. The spatial distributions indicated by the two surveys are 
mutually consistent at the angular scale of the LAT PSF except for the overall normalization. 

The first H 2 template, I^-template- 1 , is made by combining the two surveys and accounting 
for their relative intensity scales (Fig. 2b): NANTEN Woo for the area defined by the solid white 
line and that by Dame et al (2001) for the rest of the region. We refer to the 3 regions defined 
by dashed lines in Fig. 3b as “the 3 Orion regions" hereafter 9 . Wc scale the NANTEN data by 
the factor 1/1.11 to adjust the intensities to a common scale with Dame et al (2001) because the 
survey by Dame et al. (2001) has been widely used in gamma- ray analyses. 

We first start tbe analyses by setting one common D % factor for Woo in the ROI (Ha-template- 
1), or equivalently, one common Xco for the entire ROI. In the second H 2 template, Hj- template-2, 
the Wco distribution is divided into 4 regions (the 3 Orion regions and the rest of the ROI) and 
allow B x , or equivalently Xco 1 to be different in each region. We add a "dark gas” template 
(Grenier et al 2005) to H^tezn plate- 1 to make the third H 2 template, H2-tcinplate-3 (Fig. 2a and 
2c). The normalization is set free for the 2 templates. 

These spatial templates are described further in the subsections to follow. 


8.1.3. Inverse Compton template 

The inverse Compton component is known to be minor around the Orion clouds. We use the IC 
spectrum and spatial distribution given by GALPROP where the interstellar photon fields arc taken 
from Porter et al. (2008). The typical Galactic-scale IC intensity in the region is ~ 5 times smaller 
than the isotropic component described later, and their spatial and spectral distributions are si milar 
in this region. Possible local enhancement is the IC emission around the Orion Nebula (M42) 
where strong ultraviolet emission (e g. Murthy et al. 2005) and moderate infrared emission (e.g. 
Priainzano et al. 2008) exists. According to our calculation, such IC emissions are not detectable 
with the current LAT sensitivity (Orlando Ic Strong 2008). 


’The boundaries are Orion A Region I (217° > / > 212°. -23° < b < -16*), Orion A Region II (212* > £> 205°, 
-23 c < b < -16°. excluding the overlap with Onon B), and Orion B (209 6 > l > 203°, -18* < b < -13°) 



3.1.4. Point sources in the Orion region 


More than 1400 point sources are reported in the First Fermi LAT Catalog (Abdo et al 
2010b). Among them, 30 point sources are in our ROI, (l,b) = (210 ± 15°, —20 ± 15°). There 
are an additional 29 sources within 5 deg of the ROI. In the likelihood fit to be discussed later, 
the normalization is set free, energy-bin by energy-bin, for 25 high-confidence sources in the 
ROI; the indexes and normalizations are fixed to the values given in the First Fermi LAI' Cat- 
alog (Abdo et al. 2010b) for those outside of the region. There are 5 low-confidence sources 
(or candidates) overlapping with the clouds: they are 10 1FGL J0540.4-0737c, J0536.2-0607c, 
J0534.7- 0531c, J0541.9-0204c, and J0547.0r0020c Their fluxes are all low and labeled as V 
in the catalog, meaning either their flux estimates are uncertain, or they can bo artifacts resulting 
from incorrect modeling of the Galactic diffuse emission. We fit the spatial templates and analyze 
the spectra in the 3 Orion regions with and without them. The results we quote will be for the 
analyses without them: we include their possible contribution in the systematic error. 


3.1.5 . Isotropic component 

In the present analyses, the extragalactic emission and residual CR background in the data are 
not separated but treated as a single isotropic component. (Abdo et al 2009c, 2010c; Ackermann et al. 
2010; Abdo et al. 2010a). The total flux of the component at 1 GcV is ^ 25% of that associated 
with Hj when averaged over the 3 Orion regions (subtending ~ 30 msr) defined around Orion A 
and B (see Fig. 3b). 

The residual background in the Pass6 Diffuse class consists of CR-induced events raisclassified 
as gamma rays and CRs that converted in the passive material just outside of the LAT without 
leaving a signal in the anti-coincidcncc detector (Atwood et al. 2009). When averaged over many 
orbits of observations, the residual background can bo approximated as isotropic 


3.2. Fit to the Spatial Distribution 

All spatial templates described in the previous subsection were convolved with the LAT ex- 
posure and PSF. The spatial fit is made using the binned likelihood program gtlike included in 
the Fermi ScienceTools* 1 and the 4 normalizations (A,. B it Cy and Dj) in Eq. (1) are determined 
independently for the 22 onergy bins. Wo note again that ICi are fixed at the values given by 
GALPROP. Each H I annulus has a separate A,. We report only A, for the local annulus as others 
are not determined well because they lie mostly outside of our ROI. 


>0 No new sources have been added in this region in the Second Fermi LAT Catalog (Abdo nr. al. 2011) 
11 We use ScienceTools version vQrlCpO wi»h P6 VS. DIFFUSE .nulrument responw? functions. 




Our scientific interest is to study tho contributions from the gas concentrations identified as 
Orion A and B, which arc believed to be predominantly H 2 . We consider, hence, the sum of the 
H I, IC, point-source, and isotropic components as the “background” which is determined by fitting 
the observed gamma-ray distribution for each of the 22 energy bins. In the fits, we assume that H 2 - 
tcmplate-1, or the Wco distribution, represents approximately the Hj distribution. The gamma-ray 
distribution associated with the H 2 gas can be extracted less dependcntly on yet-unknown U 2 -Wco 
relation by subtracting the “background” from the observed gaxnma-ray distribution. 

We define 2 improved H 2 templates, H2-temp'.ate-2 and 3 after the initial analysis on H 2 - 
template-1. The spatial distribution is not proportional to Wco for the 2 improved templates and 
hence the “background” is different for each ^-template by a small amount. The difference is 
however negligible. 


3.2.1. Spatial Fit with Wco of One Xcq: ll 2 -templatc-l 

We use H^template- 1 as an approximation for the Ha gas distribution and fit Eq. (1) to 
determine the “background”. The energy-summed g&mma-ray distribution after subtracting the 
“background” is shown in Fig. 3a and that of the Wco-baaed model, or the product of E B x in 
Eq. (1) and H 2 -tcmplate-l, is given in Fig. 3b. The two count distributions are correlated pixcl-by- 
pixei (1 x 1 deg 2 ) in the 3 Orion regions in Fig 4a. We expect a good linear correlation between 
the two if Wco is a good tracer of H 2 . 

We note first that the correlation is fairly linear and gives a correlation coefficient 12 of 0.93. 
We then note that the correlation significantly improves if we separate the Orion clouds into the 
3 Orion regions, Orion A Region I (black solid line) and II (red dashed line), and Orion B (blue 
dotted line). The correlation coefficients for the 3 Orion regions are 0.98, 0 96, and 0.98, and the 
best-fit slopes are 0.72, 0.99, and 1.25, respectively. 

The large difference (~ 40 - 60%) in the best-fit slope suggests that the mass column density 
in Orion A and B cannot be simply derived using the same value of Xqq. We find more gamma 
rays in Orion A Region I per Wco than in Orion A Region II and Orion B, suggesting Xco is 
different in the 3 Orion regions, or that some fraction of the H 2 gas is not traced by Wco provided 
a uniform CR density. We explore these two possibilities by redefining the H 2 -template. 


3.2.2. Spatial Fit with Wco of 4 different Xqo values: H 2 - template- 2 

Based on the relation found between the spatial distributions of the gamma-ray intensity 
associated with the H 2 gas and the Wco- based model (H 2 -template-l), we make a second template. 


’’The correlation coefficient iB defined as E(x - *)(y - - f ) 3 





Ha-tcmplate-2, that will delineate tho H 2 column density more faithfully. In the template we divide 
the ROI into 4 regions, the 3 Orion regions and the rest of the ROI, and allow B, to be different in 
each region, or introduce 4 B,’s. 

The fitted results for A\ (H I) and B, (H2-template-2) in Eq. (1) are listed in Table 1 after 
combining the highest 10 energy bins into 3 bins. The gamma ray count map is shown in Fig. 3c 
is tho sum of the 4 B,’s multiplied with the corresponding components of H2-template-2. We note 
that the 3 Orion regions mix to some degree through the Fcrmt PSF. The correlation between the 
gammarray distribution associated with H 2 and the H 2 template improved as shown in Fig. 4b: the 
best- fit slopes for Orion A Region I, Region II and Orion B are 0.95, 0 94, and 1.03, respectively, 
while the correlation factors remain almost the same, 0.98, 0.99, and 0.96, respectively. 

The Xco for the 4 regions can be calculated directly as the ratio of B» to 2A t (the H 2 /H 1 
method) or by extracting the gamma-ray emission in the regions (the pionic method). The results 
from the former are given in Table 2 together with those from tho latter which will be described in 
Subsection 3.3. 


8.2.$. Spatial fit with Wc o and “ dark gas”: H 2 -template-3 

Grenier et al. (2005) found that a significant fraction of local diffuse gamma-ray emission 
observed by EGRET is not associated with either H 1 or Woo, but rather with the dust map 
traced by thermal infrared emission given by Schlegel et al. (1998). The missing gas component Is 
often referred to as the “dark gas” Other LAT observations have found gamma rays associated 
with such “dark gas” (Abdo et al. 2010c; Ackcrmann et al. 2010). We note recent measurements 
of attenuation or reddening of background stars have also detected gas concentrations not traced 
well by Wco (Dobashi et al. 2005; RowIkh h Froebrich 2009; Dobashi 2011; Ade et al. 2011). 

We make a third template, ^-template 3, that can briug out the true gas distribution associ- 
ated with the Orion clouds and enhance our understanding of the Wco-to-Ha relation by introducing 
the “dark gas". The new H 2 template consists of H?- template- 1, or Wco. and a “dark gas' spatial 
template with a normalization factor for each. 

Our “dark gas" template has bocu produced following the prescription given by Grenier et al. 
(2005) and referred to as £(B - V) r— . It is a residual map obtained by subtracting the best- 
fit linear combination of N(H 1) and Wco from the E(B - V ) map of Schlegel et al. (1998) as 
described in Ackcrmann et al. (2010). Fig. 2c shows the E{B - V) TW map around our ROI. There 
is a problem with the color temperature correction of the map by Schlegel et al. (1998) around 
the OB associations in the Orion A and B clouds, and thus E[B - V)r« value is negative In these 
points. We masked out these pixels in the E(B - V)r« map by setting the corresponding values 
to zero. 


The results for A x (H 1) and B, (2 normalizations, one for Wco and one for the “dark gas") in 



Eq. (1) are listed after combining the highest 10 energy bins into 3 bins in Table 3. The distribution 
of the gamma-ray counts associated with Hj- template-3, the sum of the counts associated with Wco 
and the “dark gas", is given in Fig. 3d. The correlation between the extracted gamma- ray counts 
and the model counts improves as shown in Fig. 4c, bringing the correlation coefficients to 0.99, 0.99, 
0.97, and 0.98, for Orion A Region I, Region II, Orion B and the sum of the 3 regions, respectively. 
The improvement in the correlation, or oquivalently in the spatial fit, comes from inclusion of 
E(B - V) T ta which has the largest contribution in the Orion A Region I seen in Fig. 2c. 


3 8. 4- Summary of the Spatial Fits 

The relative likelihoods of the spatial fits with Eq. (1) in the ROl are compared among the 3 
H 2 -templates in Fig. 5 for the 22 energy bins. The “dark gas" template (H 2 -tcmplate-3) gives the 
best fit in almost all energy bins and the 3-Aco template (H? template-2) gives the second bait 
result. The improvements relative to H 2 -template- 1 arc statistically significant. 

The residuals of the fits with the 3 templates in the ROI are given in Fig. 6. The rectangular 
boundaries of the 3 Orion regions shown in Figs. 3b, c, d are replicated in the figure. The residuals 
are significant within the Orion regions for H 2 -template-l (Fig. 6a) but not for the other 2 templates 
(Fig. 6b, c), which is consistent with the improvement we saw in Fig. 4. The difference in the 
residuals for Hj-template-2 and Hj-teroplato-S in the Orion regions is not significant relative to 
the systematic uncertainty discussed in the next subsection. We find that the large improvement 
H 2 - template- 3 has brought relative to H2-template-2 in Fig. 5 comes primarily from outside of the 
3 Orion regions, especially in the Monoceron R2 region and in the northern region adjacent to the 
Orion B: the template adds “dark gas’ in that part whereas the other templates only modify the 
3 Orion regions. 

The value of Xqo has been calculated by the H 2 /H I method by taking the ratio of B{ to 
2 Ai for the parts associated with Wco »n the ^-templates and listed in Table 2. In the pionic 
method of evaluating Xco, however, the pionic component must be extracted out of the gamma- ray 
spectrum associated with the Ha-template as will be described in Subsection 3.3. We will discuss 
the systematic errors in evaluating .Yco and possible interpretations of the results in Section 4. 


3.3. Analyses of Spectra 

The spectra associated with the II 1 and H 2 - template- 1, with the H I and H 2 -template-2, and 
with the H I and H2-template-3 arc obtained by assembling the fitted results for the respective 
templates, and B t , as shown in Figs. 7, 8, and 9, respectively. The spectra are fitted as a sum 
of the pionic and bremsstrahlung components. The gamma-ray spectra associated with the spatial 
templates (H I, inverse Compton, isotropic. And sum of Xco * Wco) plotted for the 3 Orion 
regions in Fig. 10a and b. We analyze for the gamma rays associated with the 3 H 2 -templates in 



this subsection. 


8.8.1. Fit with Gamma-ray Emission Mods Is 

The spectral template of pionic gamma rays has been calculated by convolving the gamma- 
ray inclusive cross-section for p p interaction parameterized by Kamae et al. (2006) and the CR 
proton spectrum predicted by GALPROP at the. Orion clouds. 13 The proton flux is predicted in 
the Orion clouds (R 8 8 kpc, Z = —0.14 kpc) to be — 8% smaller than that at the solar system 
( R = 8.5 kpc, Z — 0.0 kpc) where the GALPROP proton spectrum has been determined by the CR 
data taken at the Earth. The value at the Orion clouds is consistent with that determined using 
the gamma rays from the local H 1 (Abdo et aL 2009c). The good fit to the data seen Figs. 7, 8, 
and 9 supports GALPROP’9 prediction of CR spectral shape in the Orkin region and the overall 
modeling of Eq. (1). 

Bremsstrahlung emission induced by CR electrons interacting with gas is calculated in GAL- 
PROP using recent bremsstrahlung calculations (Strong it Moskalenko 1998; Strong et al. 2000, 
and references therein). The electron injection spectrum in our GALPROP calculation had been 
adjusted to reproduce, approximately, the power-law index of the electron spectrum measured by 
the Fermi LAT (Abdo ot al. 2009e). In addition, the normalization of the spectrum is adjusted to 
reproduce the LAT observed ganuna-ray flux at a low-energy band. In the spectral fits described 
below, we kept the electron-to-proton ratio, or equivalently the bremsstrahlung- to-pion ratio, fixed 
to the value given in GALPROP. When we refer to the gamma-ray emissivity per atom or molecule, 
we do not differentiate the underlying processes, but rather the sum of the bremsstrahlung and 
pionic contributions. 

The spectral fit of the H I component is reasonable for all 3 templates (x 2 = 17.7, 9.9, and 
17.1 for /dof * 14 respectively) as shown In Figs. 7a, 8a and 9a. Our pionic flux associated with 
H I is consistent with that obtained in the Fermi study on the local interstellar gas (Abdo et al. 
2009c) as overlaid in Fig. 8a. We note however that there may be a small offset between the two 
as will be discussed later. The spectra associated with molecular clouds are also fitted well by the 
3 H 2 - templates as shown below. 

The mass-to-Wco ratio, JVco. can be obtained by comparing the assumed pionic gamma-ray 
emissivity per H atom with the observed gamma-ray emissivity per VFco as shown in Figs. 7b. 


l> In GALDEF 54.77Xv»rh7S, the CR proton flux «u artificially multiplied by 1.15 to reproduce gamma-ray 
observations by Fermi The factor originates from the underestimate of gamma-ray emissivity for He and heavier 
atom* in the interstellar medium (ISM) in GALPROP. Instead of using the 1.15 correction factor, we combined the 
calculation by CtuMinr k Schaefer (1992) for contributions from CR He and heavier atoms, and the calculation by 
Mori (2009) for heavier atoms in the ISM. Hence the total gamma-ray emissivity per H atom is 1.70 times larger 
than that for p - p collisions only The difference between the total gamma-ray emissivity in the two literatures is 
~ 5%, which is taken into account in the systematic uncertainty. 




8b-d and 9b The former is calculated in the unit of MoV l a" 1 sr" 1 , and the latter is measured iu 
the unit of MeV -I 3 -, sr" 1 (2 x 10 20 cm“ 2 (K km s" 1 ) -1 )” 1 . Thus, Xco/2 of the clouds is derived 
by dividing the latter by the former. 

The results of the spectral fit for the U I component are not used to determine Xco in the 
pionic method. The fits to the spectral components shown in Figs. 7a, 8a and 9a are only to check 
overall consistency of our analyses. Their normalizations are consistent within the uncertainty in 
the H I column density discussed in Section 4. 


3.3.2. Spectra obtained with H 2 -template- 1 

The fitted spectra are plotted as stuns of pionic and bremsstrahlung emissions in Fig. 7a, b 
for the H I spatial template and the ^-template- 1 (Orion A Region I, II, and Orion B combined), 
giving x 2 /dof of 17.7/14 and 20.2/14, respectively. 

We give the Xqo value obtained from the fitted pionic spectra in Table 2. Since the fit is 
substantially poorer than those for H2-template-2 and 3 (see Fig. 5), the value should be taken just 
as a reference value. For this reason we do not quote systematic errors in the table. 


3.3.3. Spectra obtained with H 2 template-2 

The fitted spectra are plotted as sums of pionic and bremsstr&hlung emissions in Fig. 8b, c, d 
for Orion A Region I, II and Orion B, giving x 2 /doi of 14.0/14. 18.5/14, and 10.6/14, respectively. 
The Xco values obtained for the 4 regions from the fitted pionic spectra are given in Table 2. 

The coefficient Xco Is significantly higher for Orion A Region I than for other regious, con- 
sistent with the slopes obtained in Subsection 3.2 in the pixel- by- pixel correlation study. This also 
can be seen in the Xco obtained with the H 2 /H 1 method. 

We note that the fraction of the H I component in the gamma-ray spectrum integrated in the 
3 Orion regions is comparable to that associated with Wc o (see Fig. 10a). This is because the solid 
angle subtended by the Orion molecular clouds is a small fraction of our 3 Orion regions in solid 
angle and the overall mass of atomic gas is greater. 


3.3.4 • Spectra obtained with H 2 -tcmplate-3 

The fitted spectra integrated over the Wfoo and “dark gas" components are shown in Fig. 9b, 
c. We give Xco for the ROI from the fitted pionic spectrum in Tkble 2. 

The Xco obtained in fits with the Wco can be compared with these obtained in similar 



analyses Including the “dark gas" template: 2.0 x 10 20 (in the local arm), 1.9 x 10 20 (the Perseus 
arm) and 0.87 x 10 20 (the Gould Belt) in the same unit as above ( Ackermann et al. 2010; Abdo et al. 
2010c). 

The spectrum aasociated with the “dark gas" component la similar in shape to that associated 
with IVqo hut about half as intense (Fig. 10b). The 2 spectral cuergy densities (SEDs) become 
comparable in Orion A Region I aa seen in Fig. 10c. The “dark gas” dominates over Wqq in the 
pixels near the high longitude end of Orion A and eventually Wco diminishes in the pixels beyond 
them towards higher longitude. 

Our Xco measurements given in Table 2 can be compared with those determined using 
the gamma-ray flux from the Orion-Monoceros complex measured with EGRET: (1.35 ± 0.15) x 
10 20 cm' 2 (K km s -1 )" 1 (Digcl et al. 1999). We note there were no Galactic CR propagation mod- 
els such as GALPROP nor CR measurements as precise as are available now: Xco determined 
by the Hj/H I method and it compares well with the single Xco value of 1.36 ±0.02 obtained with 
the H2- template-!. 


3.4. Total masses of Orion A and B 

The distance from the Sun to the Orion nebula (M42) inside the Orion A has recently been 
measured by parallax to be 389tJJ pc (Sandstrom et al. 2007), 414 ± 7 pc (Mentea et al. 2007), 
437 ± 19 pc (Hirota ct ah 2007), and 419 m 6 pc (Kim et al. 2008). We adopted 400 pc as the 
distance to the Orion A and B clouds and used the total pionic gamma-ray fluxes obtained above 
to get the total masses of Orion A and B outside 14 . 

Mass estimation using Ha tcmplate-2: 

M A = (74.5 ± 1.3) x 10 3 M 400 
M b = (33.5 ±0.7) x 10 3 M 400 

where 



and d is the distance to the clouds. We will discuss the systematic uncertainties in the next section. 

Mass estimation using H2-template-3: Addition of the “dark gas" changes the estimation of 
the Orion A and B masses by about 10%. 

A*A,Wco * (55. 1 ± 0.8) X ID 3 oo 


M W e not* that the apatlal extent of Orion U defined here is aiguificantly different from that used in Wilson et nl 
(2005) because we are unable to separate Orion B from the complex doud structure* behind due to the broad PSK 
of the LAT. 



M AfD *r k = (27.6 ± 0.7) x 10 3 M 400 
Mq = (36.0 ± 0.5) x 10 3 M 400 . 

The total mass of Orion A (== M\,Wco + ^A.Dvk) la (82.7 ± 1.1) x 10 3 M 400 The Orion A mass has 
been estimated by Wilson et al. (2005), assuming Xco = l-8x IO 20 cm _2 (K km a" 1 ) -1 (Dame ot al. 
2001), to be M\ = 91.7 x 10 3 M 400 . The maas has been estimated separately for Orion A Regions 1, 
2, 3, and NGC 2149 In Wilson et a!. (2005). Our Orion A (Region 1 and II) includes their Regions 
1, 2, and 3 but overlaps only partially with NGC 2149. Considering the breadth of the PSF and the 
limited statistics of the data, we could not determine how much of NGC 2149 overlaps our Orion 
Region I. If we assume about one half of NGC 2149 is in our Orion Region 1 and the systematic 
error introduced by this ambiguity is half of the NGC 2149 mass estimated by Wilson et al. (2005), 
the Orion A mass to be compared becomes M A = (86.3 ± 5.4) x 10 3 M 4 on- The Orion B region is 
more complex and such a comparison Is very difficult. 


4. Discussion 

Although the Orion clouds lie away from the Galactic Plane and subtend relatively small solid 
angle, many Galactic and extragalactic sources contribute to the ROI through the large PSF of the 
Fermi - L AT. 

We have analyzed the observed data to extract the intensity associated with the molecular 
clouds, the 3 Orion regions in particular, by using the 3 Hj-tcmplates made from Woo on the 3 
different assumptions for each of the 22 energy bins. The ratio of the normalization factors for 
H I and Ha, AJ'2B it gives the conversion factor of Woo to the maas column density, Xco (the 
Ha/H I method). For this, the H I mass column density must be well understood from the radiative 
transfer of the H I line and the CR spectrum must be constant in the ROI. 

In the second method (the pionic method), Xco ia determined by comparing the observed 
pionic gamma-ray intensities with those expected from the CR spectrum at the Orion clouds and the 
pionic gamma-ray production cross- section. For this, we have to know the absolute CR spectrum 
and flux, the instrument response function (IRF), and the pionic gamma-ray production cross- 
section, in particular the pionic gamma-ray contribution from metals in CR and ISM. 

In the subsections to follow, we evaluate uncertainties and possible systematic errors in the 
nnalyses, especially in evaluating Xco > n the 3 Orion regions. We then summarize the results 
obtained in this paper and present possible interpretations thereon. 


4.1. Possible Systematic Errors in the Analyses 


Systematic errors that affect tho correlation measurements between gamma-ray intensities and 
Wco are discussed in two categories: the first one applies commonly to the 3 Orion regions and 



th* second affects the relation differently in the 3 regions. 


4.1.1, CR intensity at the Orion clouds 

Uncertainty in the fluxes and spectra of CRs, in particular those of protons, can affect in both 
categories. The Galactic CR protons that produce pious in our energy range remain in our Galaxy 
longer (~ 5 x 10 7 yrs) than electrons (~ 7 x 10 6 yrs) (Lee ct al. 2011) and their flux variation within 
the Galaxy is believed to be predicted well by GALPROP. We note that the CR source distribution, 
the Galaxy size, and the CR diffusion coefficient are the important inputs to GALPROP. Using 
the CR spectrum measured at the Earth, we have calculated the CR spectrum in the Orion region 
for the 2 choices of the CR source distributions and the 3 choices of Galactic halo heights (2, 4, 
and 10 kpc) used in a GALPROP-bascd study by loc et al. (2011). The CR spectrum does not 
change more than ~ 2% from the value used here as long as it is constrained to the measurements 
at the Earth and to reproduce the Galactic diffuse gamma- ray intensities measured by the Ferm » 
LAT (see Lee ct al. (2011)). We also note that the gamma-ray spectrum from the local H 1 (typical 
distance < 1 kpc) is consistent with the CR proton flux being within ~ 10% of that at the Earth 
(Abdo et al. 2009c). 

CRs could be accelerated in the clouds and/or prevented from penetrating into their cores by 
embodded magnetic field. Wc first note that there are no strong non thermal X-ray source uor radio 
SNR found In the clouds (Feigelson et al. 2002, and references therein). Therefore no appreciable 
CR acceleration is likely to be taking place In the Onon clouds. The good linear correlation between 
Wqo and gamma- ray intensity seen in all 3 Orion regions (Fig. 4a) confirms that the CRs effective 
In producing pious (kinetic energy > 1 GeV) are penetrating well inside the higher-density parts 
of the clouds. 

Based on these observations we assume that the CR flux in the Orion region is 8% lowci than 
that at the Earth with possible systematic error of :rl0% due mostly to disagreement among recent 
CR measurements at the Earth and solar demodulation uncertainties. 

Uncertainty in the CR flux at the Orion clouds contributes directly to the systematic error in 
the pionic method but indirectly in the Ha/H I method. In the former, the absolute CR intensity 
is assumed to be known while the CR intensity is assumed to be the same in the local II I region 
and the molecular clouds in the latter 


4- 1.2. Uncertainty in the instrument response functions 

The uncertainty in the absolute calibration of the LAT effective area can also introduce error 
of the first kind. The effective areas were derived based on Monte Carlo studies of the LAT, checked 
against beam tests at accelerators (Abdo et al. 2009a; Atwood ct al. 2009). Comparisons between 



flight data and Monte Carlo studies have been made to quantify the systematic uncertainty in 
the effective area (Abdo et al. 2009b). At present, we estimate this systematic error to be 10% at 
100 MeV, 5% at 500 McV and 20% at 20 GeV. 

The systematic error in the absolute energy scale has been estimated as -1-5/ — 10% (Abdo et al. 
2009c). We have refitted Xco after artificially shifting the energy scale by +5% and by - 10%: the 
number of pionic gamma rays changes less than +1/ — 8% for all 3 Orion regions with all 3 H2 
templates. Wo include this possible error due to the uncertainty in the energy calibration when 
assessing the overall systematic error. 

The pionic method is affected directly by the uncertainty in the instrument response function 
while the H2/H I method is insensitive because it affects the denominator and numerator similarly. 


4.1.3. Uncertainty in the spin temperature of H I 

In converting the observed 21 cm line emissiou intensity (Kalberla et al 2005) to the H 1 
column density, Tg was assumed to be 125 K. The range of Tg measured in the local H 1 gas varies 
broadly between 90 K and 400 K (e g., Mohan et al 2004a,b, and references therein) while we 
have assumed a likely range for our ROI to be between 90 K and 250 K. 

We refitted the Fermi data in the ROI with these two extreme Tg values with ^-template 
2 and 3. We then calculated Xco by dividing B % by 2A* in Eq. (1), or by extracting the pion 
component in the spectra. The deviations of Xco from those obtained with Tg of 125 K are taken 
into account in the systematic errors given in Table 2. The large systematic errors for Xco on 
B t /2At (Column 3) enter via 2A, which depends on the absolute calibration of the H I gas density 
or Tg in the local H I. The pionic method uses the product of the CR intensity and pp -► 7 cross- 
section in place of 2 A, and is less directly affected by the uncertainty in H I gas density or Tg of 
the local H I, although the uncertainties can have a small indirect effect through the overall spatial 
fitting. This effect is much smaller than the overall systematic error and negligible. Wc note that 
there is some discrepancy between the gamma-ray spectra associated with H I in the ROI and the 
local H I (Abdo et al. 2009c) as seen in Fig. 8. 


4.I.4. Effect of overlapping point source candidates 

We have not included the 5 sources overlapping with the Orion clouds (Sec. 3.1.4) because 
they are all classified as “potentially confused with interstellar diffuse emission or perhaps spurious" 
(Abdo et al. 2010b). To investigate their potential contribution we repeated the analysis including 
these sources with the fluxes and spectra listed in the First Fermi LAT Catalog. The fit with the 
pionic method gives the following Xco ia unit of cm -2 (K km s" 1 )” 1 : (2.29 ±0.05) x 10 20 for Orion 
Region I; (1.16 ±0.05) x 10 20 for Orion Region II; and (1.24 ±0.04) x 10 20 for Orion B. They are 2%, 



19%, and 8% less than those obtained without these point source candidates. In the present study, 
we assume they are artifacts and add +0/ — 2, +0/ - 19, and +0/ - 8% to the overall systematic 
error in the 3 regions. 


4-1.5. Overall error 

For the Hj/H I method, the uncertainty in the H I mass density 20%) due mostly to the 
uncertainty in 7’s dominates the systematic error . Other contributions include the overlapping V 
sources (-r0/ - 2, -t-0/ — 19, and +0/ - 8%) and variation in the CR intensity within ~ 1 kpc or 
between H I and the molecular clouds (±5%), making the total systematic errors for the 3 Orion 
regions to +25/ — 28, +25/ - 44, and —25/ - 33% as given in Column 3 of Thble 2. 

For the pionic method, the overall systematic error in determining Xco comes from the un- 
certainty in the IRF including that, due to the energy calibration uncertainty (±10%), unknown 
contributions of the overlapping sources (+0/- 2, +0/- 19, and +0/ - 8%), uncertainty in the CR 
intensity (±10%), uncertainty in the pp pion production cross- section (±5%), and uncertainty in 
the contribution from heavier nuclei (±5%). We conservatively quote the linear sum of these com 
binations as the possible systematic error for the 3 Orion regions, which are +30/ - 32, -t 30/ - 49, 
and +30/ — 38%, as given in Column 5 of Table 2. 

The systematic errors that can affect Xco differently in the 3 Orion regions are variation in 
the CR intensity within ~ 1 kpc (r:5%) and the overlapping sources. The overall error of this kind 
is conservatively estimated to be the linear sum of the two, +5/ 7, +5/ - 24, and +5/ - 13%. 


4.2. Gamma-ray intensity and E(J — H) 

The line-of-sight visual attenuation, 4v, are often used as a gas mass tracer in theory- based 
studies of the CO fraction in all molecules including carbon and hydrogen (e.g., Burgh et al. 2010; 
Wolfire ct al. 2010; Glover et al. 2010, and references therein). To calibrate crudely our mass col- 
umn density with A\ used in these theory-based analyses, wc have related the gamma-ray counts 
on the horizontal axes of Fig. 4 and E ( J - H) in the 3 Orion regions measured by Doba&hi (2011). 
We note that the atomic and molecular components are assumed to be contained within a fixed 
length (e.g. 20pc) along the line-of-sight in the theory- based analyses while the components are 
measured as column densities integrated over unknown lengths along the line-of-sight in observa- 
tions. Moreover E(J H) is known to trace the Ha gas but also pick up some H 1 gas through dust 
mixed with it. Hence the cross calibration works at best crudely and only in the regions of clouds 
where the H 2 longitudinal distribution is well confined and the Ha volume density dominates over 
that of H I. Despite these uncertainties, it is important that our measurements be compared with 
theory-based analyses. 



We found good linear relations for the pixels with high gamma-ray counts (> 300 per deg 2 ) 
in all 3 Orion regions and could correlate the gamma- ray count scale on the horizontal axes of 
Fig. 4 to E{J — H) assuming Ay * Rv-ejh * E{J — H). The Rv-ejh has been determined 
observationally and its value ranges between 7.8 (Dobashi 2011) to 10.9 (Cardclli et al. 1989). The 
highest point in our count map is ~ 700 per pixel in Orion A Region II where Ila concentration 
is highest and the corresponding value of A\ is ~ 5 when averaged over 1 x 1 deg 2 pixels for an 
assumed value of Ry ejh = 7.8. So A\ = 5 on the horizontal axes of Figs. 5 and 6 in Glover et al. 
(2010) corresponds crudely to ~ 700 counts per pixel assuming H 2 is well confined (e.g., to ~ 20 pc) 
along the line- of- sight. 


4.3. Summary of the Results 

The results obtained in the present work are significant beyond the estimated systematic errors. 

They are: 

1 Linearity holds between mass density associated with the Orion clouds and Wgq\ Ab discussed 
in Subsection 3.2 and shown in Fig. 4, our results suggest that CRs penetrate to all translucent 
part of the clouds. Possible shielding of CRs discussed in Aharonian (2001) does not apply 
to most parts of the Orion clouds. 

2. The Xco factors calculated with tho piouic method and with the Hj/H 1 method differ by 
~ 15% but agree within the estimated systematic error (Table 2). The difference can be 
explained by uncertainties in the column densities of H 1 and calculation of gamma-ray emis- 
sivity per H I atom. 

3. The Xco factor obtained with the H2-templaxe-2 is found to be larger by ** 40 - 60% in 
Orion A Region I than Orion A Region II and Orion B for the two methods. The difference 
is much larger than the systematic error that can affect the Xco factor differently in the 3 
Orion regions (Thble 2). 

4. In the “dark gas” scenario, the added “dark gas" accounts for the majority of the gas not 
traced by Wco- One Xco factor can then describe the Wco-traced H 2 distribution in the 
ROL 


4.4. Interpretation of our results on Xco 

Historically the relation between A’(H 2 ) to Wco has been considered to depend on the envi- 
ronment around the molecular cloud. The environmental factors discussed in the literature are: 

Metaliicity: This possibility has been discussed in the literature since the late 1980'a (e.g., Elmegreen 
1989; Bolatto et al. 1999). According to an empirical formula proposed to relate Xco to 



[0/H] (Wilson 1995; Arimoto et al. 1996), the metallicity must be ^ 2 times higher in Orion 
A Region I to account for the observed difference in Xco between Region I and II, which is 
unlikely according to Galactic-scale measurements (e.g. Esteban ct al. 2005). We note that 
metallicity is generally considered to be an important environmental factor influencing the 
Ha-to-H I ratio. 

Overlapping H I clumps: Compact H I clouds with angular diameters of 1 — 2dcg have been 
found in various Galactic locations (e g. Braun k Strom 1986; Kavars ct al. 2003; Lee et al. 
2008). A new reanalysis of the LAB H I survey shows no such concentration detected at the 
sensitivity level of the present study overlapping with the Orion A and B (Kalbcrla ct al. 
2010 ). 

Low density H2 not traced well by Wqo’ Existence of diffuse H2 gas not traced well by Wco 
has been discussed in the literature cited in Section I and Fermi analyses are bringing the 
discussion to a quantitative level (Ackermann et al. 2010). We refer to the following recent 
works on the Ha and CO fractions and try to interpret our results: 

• Burgli ct al. (2010) have studied the fractions based on Hubble Space Telescope, observa- 
tions and characterized the Xco dependency on jV(Ha). 

• Wolfire et al. (2010) have studied chemical composition of a model cloud theoretically 
and found that CO becomes depleted because of photodissociation in the periphery 
where the gas density decreases. 

• Glover ct al. (2010) have studied the time-dependent H2 and CO fractions in clouds 
through computer simulations and found Xco increases sharply where jV(Ha) decreases 
for Ay/ < 3.5. 

All of the above studies predict that the CO/(total C) fraction drops as the H2 column density 
decreases, as toward the periphery of Orion A and B. However the Wco*fr>Ha relation and 
the abundance of Ha-without-CO gas may be more complicated Fbr example, Ikcda et al. 
(2002) found that Af(C I)/Af(CO) increases to high values along all of the peripheries whereas 
we find Region I of Orion A to be more abundant in CO-depletcd gas than Region 11. The 
prediction that Xco increases sharply in regions Ay/ < 3.5 by Glover et al. (2010) Is consistent 
with our finding that the “dark gas" is concentrated in the high- longitude end of Orion A 
where Wco becomes low. 


5. Conclusion 

We have reported on the first 21 months’ observations of Orion A and B with the Fermi 
Gamma-my Space Telescope in the energy band between ~ 178 MeV and ~ 100 GeV We have mea- 
sured the mass column density distribution within the clouds at the angular scale of the instrument 
PSF using the pp -► 7 production cross-section accurately calibrated at accelerators as well as using 



the gamma-ray emissivity of the local H 1 gas. We found with the pionic method that a linear rela- 
tion holds between mass density and Wco with Xco — 2.34, 1.43, 1.35 x 10 20 cm _2 (K km a -1 ) -1 
with a systematic uncertainty of +5/ — 7, +5/ - 24, and +5/ - 13% (relative in the 3 regions), and 
+30/ — 32, +30/ — 49, and +30/ — 33% (absolute) for Orion A Region I, Region II, and Orion B, 
respectively. Those values are consistent with the Xco values determined with the more traditional 
H 2 /H I method (Xco = 1.97, 1.20, 1.14 xlO 20 cin -2 (K km s” 1 )” 1 ) within our overall systematic 
error. This implies that Galactic CRa arc penetrating into most parts of the clouds. The analyses 
also included the “dark gas” (Grenier et nJ. 2005) not traced by CO or H I. We found that the 
gamma-ray flux associated with the “dark gas” spatial template exceeds that associated with the 
Woo template in Orion A Region I. The situation is reversed in Region II and in Orion B. This is 
generally consistent with the fit finding t higher Xco value for Orion Region I in the absence of 
the dark-gas template. 

We have interpreted the increase in Xco *nd “dark gas” fraction in Orion A Region I in the 
light of recent studies of the relation between the H 2 and CO fractions by Burgh et al. (2010); 
Wolfire et al. (2010); Glover et al. (2010). Xco ia expected to increase rapidly as the gas oolumn 
density decreases to Ay ~ 3.5 or less (Glover ct al. 2010). The mass column density wc have 
measured in Region I corresponds to Ay < 4, dose to the predicted threshold for onset of the non- 
linearity predicted between Wco and JV(Ha). The mass column density drops further (Ay < 2) 
toward the high Galactic longitude end of the Orion A where the gas becomes “dark” to Wco. 
consistent with the predicted non- linear relation. 

The Fermi-LAT collaboration is continuing to reduce uncertainty in the IRP, identify extended 
gamma-ray sources, and improve the modeling of the Galactic-scale diffuse gamma-ray emission. 
We expect the systematic uncertainties quoted in subsection 4.1 to be reduced significantly through 
these efforts. The systematic uncertainty in the OR spectra and the H 1 mass density also will be 
reduced when the data from new experiments and surveys become available. The present analyses 
can then be updated to a higher precision and tho relation among Wco and the gas mass density 
characterized further for various molecular clouds in the Galaxy. 
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T^ble 2. Xco obtained on H 2 -templato-l, 2, and 3 


Region 

Xco a on D/2A 

SyB. error 6 (%) 

Xco a on pion 

Sys. error* (%) 

H 2 - template- 1 

Entire ROI 

1.36 ±0.02**, 

NA 

1 .63 + 0.02** 

NA 

H2-tempiate-2 

Orion A Region I 
Orion A Region IT 
Orion B 
Elsewhere 

1.97 ±0.05** 
1.20 ±0.03*^ 
1.14 ±0.03** 
1.43 ±0.04** 

+25/ - 28 
+25/ -44 
+25/ -33 
NA C 

2.34 ± 0.05** 
1.43 + 0.04** 

1.35 ± 0.03** 
1.69 ±0.04«* 

-r30/ - 32 
+30/ -49 
+30/ -38 
NA d 

H2-template-3 

Entire ROI 

1.21 ± 0.02** 

+25/ - 37%« 

1.32 ± 0.02,i* 

+30/-4tr 


•In unit of 10 20 cm~ 2 (K km s -1 ) *. 

b The systematic error is discussed in Subsection 4.1: it comes from a combination of 
uncertainties in the H I spin temperature and in the fitting process. The systematic errors 
which may apply differently to the 3 Orion regions are +5/ - 8, +5/- 24, and +5/ — 13%, 
respectively. 

c Tbe systematic error is discussed in Subsection 4.1. The systematic errors arc the same 
as b. 

c We have not attempted to estimate systematic, error outside of the Orion regions in 
this study. 


r The average of the systematic errors estimated for the 3 Orion regions. 









Table 3. Gamma-ray emissivity fitted with H2-tcmplate-3 


Energy R»n«o (MeV) Err.urvity per H •tom* Emlmiviey per Wco h Eraiseivity per E(B - V )«»' 

178-237 (4.61 £ 0.08) x 10 * (6.56±0 31) x lO" 38 (1.00 £ 010) X IQ-* 7 

237 - 316 (2 99 £ 0.08} X 10"* (3.39 £ 0.27) X lO" 38 (5.60 £ 0.83) x 10" 3 " 

316-422 (I 68 ±0.07) x 10“" (2.03 * O il) X lO"** (3.96 £ 0.36) x 10"” 

422 — 662 (1.02 £ 0.08) X lO" 30 (1 13 ±0.05) x I0" w (2 07 £ 0.16) X 10 " 

502 - 750 (5.39 £ 0 08) X ID" 30 (6.51 ±0.21) x 10 (1.37 £ 0.08) X lO" 38 

750 - 1000 (2.97 x 0.00) X 10"* (3.54 £ ai6) X 10"*° (6 57 £ 0.60) x 10" 18 

1000- 1334 (1 58 J: 0.05) x 10- 50 (1 Ml £ 0 09) X lO" 50 (3 57 £ 0.35) x 10 39 

1334-1778 (8.00 £ 1.02) x lO' 31 (9.37 £ 0.43) x 10" 31 (1 .86 £ 0.16) x 10"" 

1778 - 2371 (3.64 £ 0.25) x 10“" (5.00 £ 0.31) x 10 ^ (7.45 £ 1.19) X 10" 30 

2371 - 3162 (1.51 £ 0.14) X 10 “ (2.19 £ 0 .1T) x 10“ 31 (4 82 £ 0.68) x lO” 30 

3162 - 4217 (6 .56 £ 0.80) X 10 “ (1.06 £0.10) X 10 -•* (2.18 £ 0.40) X lO" 80 

4217 - 5623 (3 82 £ 1.71) X 10" 33 (4.31 £ 0.49) x 10 (6 50 £ 2.34) x 10~ sl 

5623- 10000 (1.06 £0.14) x 10 -33 (1-07 £ 016) X 10" 33 (2.20 £ 0 68) x 10" 31 

10000 - 23714 (1^5 £ 0.15) x 10-” (1 68 £ 0.24) x 10" 33 (1.72 £ 0.92) X 10"” 

23714 -100000 (4.62 £6 52) X 10-“ (9 55 * 3.50) x 10"“ (1 46 £ 1-22) X 10' 33 


Note. — Error* ere statistic*! only 
•MeV-’a- 1 *- 1 per M atom 
‘•MeV- J a , «r- l (2 x lO^cm" 3 ^ km s" 5 )" 1 ) 1 
‘MeV I «- 1 «r- 1 (2 x lO^rai*)- 1 
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Fig. 2. — (a) JV(H I) template summed over the line-of-sight velocity. The pixel sire is 0.5° x 0.5°, 
The dashed lines show the boundary of the ROI. (b) Wco template used in ^-template- 1 and 
Hj-template-2. We used NANTEN data (Fukui et al. 2011) in the area bounded by the solid lines 
and those by Dame et al. (2001) elsewhere. Pixel resolution is 0.125° x 0.125°. (c) E(B - V) rm 
template used in H 2 -template-3. Pixel resolution is 0.5° x 0.5°. 
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Fig. 3. — (a) Gamma-ray count distribution in the ROI after subtracting the fitted “background" 
that is the sum of the H I, IC, point-source, and isotropic components, (b) The fitted model 
map obtained by assuming one common Xco for the ROI (^-template- 1). Dashed lines define 
the boundaries of the 3 Orion regions, Orion A Region I, Region II and Orion B. (c) Samo as 
(b) but obtained by assuming 4 different Xco for Orion A Region 1, Region II, Orion B, and 
elsewhere (Ha-template-2). (d) Same as (b) but obtained by adding E{B V r ) ra to Ha- template- 1 
'Ha- tcmplate-3). 
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Fig. 4.—* (a) Correlation between the gainma-ray count distribution shown in Fig. 3& (the hor 
izontal axis) and that fitted with Ha-template-1 in Fig. 3b (the vertical axis) for all pixels in the 
3 Orion regions. Points represent pixels in Orion A Region I (black circles), Region II (red tri- 
angle), and Orion B (blue squares) with fitted lines black, red, and blue, respectively. Error bars 
represent statistical errors in counts in pixels. Same after replacing the vertical axis for that fitted 
with Ha-template-2 (b) and for that fitted with Ha-template-3 (the sum of IVco and E(B - V) im 
components) (c). 




Fig. 5. — Difference in log (Likelihood) between the spatial fit using H2-template-3 (dotted line) 
and cither that with ^-template- 1 (dashed line) or that with Ha-template 2 (solid line) In the ROI 
for the 22 energy bins. Note that the lines are drawn between the data points only to guide the 
eye. 



Fig. 6. — Residue in the energy-suznmod gammarray counts of the spatial fit with Ha-template-1 
(a), Hj-tcmplate-2 (b), and H2-t«nplate-3 (c), binned in 1 x 1 deg 2 pixels. The black dotted line® 
show the boundaries of the 3 regions, Orion A Region I, II, and Orion B. 
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Fig. 7. — Spectral energy densities (SED) associated with local H I (7s = 125 K assumed) (a) and 
that associated with fy-tcmplate-l (b). The lines are: total (solid), brems&trahlung (dotted) and 
pion decay (dashed). The CR spectral shape and electron- to-proton ratio at the Orion clouds were 
fixed to those used by GALPROP. The vertical axes are normalized to the column density of H I in 
unit of 1cm ' 1 for (a) and to 2 x Xqo in unit of 10 w cm _3 (K km s -1 ) -1 for (b). The energy bins 
between No. 13 and No.22 are combined to wider energy bins. Vertical bars represent statistical 
errors. Note that the spectral fit to H I is not used in evaluating X<x>. 
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Fig. 8. — Spectral energy density (SED) associated with local H I [T& = 125 K assumed) (a), Orion 
A Region I (b), Region II (c), and Orion B (d) for the fit with H'2-template-2. The corresponding 
SED obtained for the local H I (Abdo et a! 2009c) is shown by gray squares in (a). The assumption 
about the CR, the line legends, and the vertical axis units are the same as in Fig. 7 
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Fig. 9. SED associated with local H 1 (7s ■* 125 K assumed) (a), that associated with Woo 
(b), and that associated with E(B — V)^ (c) obtained with H^-templatc 3. The line legends and 
vertical axis units are the same as in Fig. 7. 
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Fig. 10. — Gamma-ray spectra spatially associated with two H 2 templates in the 3 Orion regions 
marked in Fig 3b: (a) the sum of the 3 regions obtained with H2-tcroplate-2; (b) the sura of the 3 
regions with H2-teinpIate-3; (c) Orion A Region I obtained with H2-template-3. Black circles show 
the isotropic component, red squares H 1, green upward triangles CO, and purple dashed hoe the 
inverse Compton. Blue downward triangles in (b) and (c) represent the spectra associated with 
E(B - V)„. 




